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We propose a one-loop induced radiative neutrino mass model with anomaly free flavour 
dependent gauge symmetry: /i minus r symmetry t/(l) /i _ T . A neutrino mass matrix satisfy¬ 
ing current experimental data can be obtained by introducing a weak isospin singlet scalar 
boson that breaks f7(l) M _ r symmetry, an inert doublet scalar field, and three right-handed 
neutrinos in addition to the fields in the standard model. We find that a characteristic 
structure appears in the neutrino mass matrix: two-zero texture form which predicts three 
non-zero neutrino masses and three non-zero CP-phases from five well measured experi¬ 
mental inputs of two squared mass differences and three mixing angles. Furthermore, it is 
clarified that only the inverted mass hierarchy is allowed in our model. In a favored param¬ 
eter set from the neutrino sector, the discrepancy in the muon anomalous magnetic moment 
between the experimental data and the the standard model prediction can be explained by 
the additional neutral gauge boson loop contribution with mass of order 100 MeV and new 
gauge coupling of order 10 -3 . 
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I. INTRODUCTION 


Radiative neutrino mass models are one of the most promising scenarios at TeV scale physics 
to explain tiny neutrino masses. The original model based on the idea of radiative generation 
of neutrino masses is known as the Zee model |l] proposed in early 80’s, where neutrino masses 
are generated at the one-loop level. After the Zee model, the Zee-Babu model [sj] has also been 
proposed, where neutrino masses are explained at the two-loop level. In 2000’s, radiative neutrino 
mass models have been extended so as to include a dark matter (DM) candidate by introducing 
an unbroken symmetry such as a discrete Z 2 symmetry known as; e.g., the models by Krauss- 
Nasri-Trodden P] and by Ma 4, jj. After these models appeared, various kinds of extensions have 
been considered in the scenario based on the raditive neutrino mass generation such as models 
with the supersymmetry [f|, the B-L symmetry Q], flavour symmetries Q| and the conformal 
symmetry 9]. Furthermore, in Refs. 10], a complex SU(2)l triplet scalar field is introduced, in 
which the collider phenomenology can be rich because of the existence of doubly-charged scalar 
bosons. Loop induced Dirac type neutrino masses have been proposed in Ref. [11]. In models 
proposed in Refs, jisj], charged lepton masses are also introduced at quantum levels in addition to 

3ove mentioned models, a number of radiative neutrino mass 


13Ml5l] up to now, and they have been classified into some groups 


neutrino masses. In addition to the a 
models have been constructed 
in Refs. flfi| . 

On the other hand, Abelian gauged U{ 1) symmetries are well compatible with such radiative 
models. It has been known that there are four different anomaly free and flavour dependent types 
of U( 1) symmetries in the leptonic sector; namely, L e — LL e — L T , and L„ — L T , where Li 


denotes the lepton number with the flavour i. Especially in the case of — L T 


17H26I]. constraints 


on the mass of additional neutral gauge boson Z' and the new gauge coupling constant from the 
LEP experiment are very weak, because the Z' boson does not couple directly to the electron. 
We thus can consider a light Z' boson scenario, by which the discrepancy in the muon anomalous 
magnetic moment between current data and the prediction in the standard model (SM) [l?J can be 
explained with the mass of Z' to be 0(100) MeV and the U( l) M _ r gauge c oupling to be 0(1CT 3 ). 


The positron anomaly reported by AMS-02 


can be probed at the 14 TeV run of the LHC 


2711 c ould be explained [13,12a]. Such a light Z' boson 


251 | through multi-lepton signals. 


In our paper, we combine a radiative neutrino mass model at one-loop level and the gauged 
U(l)fj,- T symmetry to get neutrino masses, mixings, and dark matter candidates. We find that 
a predictive two-zero texture form of a neutrino mass matrix can be obtained corresponding to 
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TABLE I: The charge assignments of leptons and scalars under SU(2)l x U(1)y and Z 2 symmetry. The 
index i(= e,g,r) denotes the lepton flavour. 
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TABLE II: The charge assignments under the gauged f7(l) p _ r symmetry. Fields which are not displayed in 
this table are neutral under U(l)fj,- T - 


“Pattern C” in Ref. [29(]. In this texture, we only need five experimental inputs to determine all 
the neutrino parameters. We can choose the most accurately measured ones: two squared mass 
differences and three mixing angles. It turns out that only the inverted mass hierarchy is allowed in 
our texture. Non-vanishing one Dirac and two Majorana CP-phases, and non-zero three neutrino 
mass eigenvalues are predicted. 

This paper is organized as follows. In Sec. II, we define our model, and give mass formulae for 
scalar bosons. In Sec. Ill, we calculate the mass matrices for the lepton sector; charged leptons, 
right-handed neutrinos and left-handed neutrinos. The detailed analysis for the two-zero texture 
form of neutrino mass matrix is also discussed. In Sec. IV, we discuss new contributions to the 
muon g — 2 and lepton flavour violation in our model. Conclusions and discussions are given in 
Sec. V. 


II. THE MODEL 


We consider a model in the framework of the gauge symmetry of SU(2)l x U{1)y x U( l) M _ r 
with an unbroken discrete Z 2 symmetry. The particle content in our model is listed in Table Q] 
The charge assignment for the U( l)^_ r symmetry is separately shown in Table HT1 

Our model is an extension of the model proposed by Ma [ 4 ], where neutrino masses are generated 
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at the one-loop level. In the Ma model, three right-handed neutrinos and an inert scalar doublet 
field are added to the standard model (SM). We introduce only one additional SU{2)l singlet 
scalar field S with the even parity under Z 2 to the Ma model. The vacuum expectation value 
(VEV) of S breaks the U( l) /i - T symmetry. 

The mass terms for right-handed neutrinos N l R and the relevant Yukawa interactions are given 
by 

-£ y = \M ee NfN% + Im^N^N^ + N^N^) + h.c. 

+ y e L e L ®e R + + y T L T L $T R + h.c. 

+ K + W n r) s * + KAW n r + W n r) s + h - c - 

+ fjj[{ia 2 )rfN* R + /„!£(« T2)v*K + N T R + h.c. (II.l) 


The scalar sector of our model is composed of a singlet ( S) and two doublets, one active (4>) and 
one inert ( 77 ). The most general scalar potential is given by 

V = ^\ 2 + ^ + y? s \S\ 2 

+ ^ll^l 4 + ^A 2 |7 ?| 4 + A 3 |4>| 2 |t7| 2 + A 4 |T t ? 7| 2 + ^A 5 [( 4 > t ? 7) 2 + h.c.] 

+ ^sI'S'l 4 + -^5<I'l'S'hl l ^| 2 + / ^5r;l‘S , | 2 | 7 7| 2 ’ (H-2) 


where all the parameters can be taken to be real without any loss of generality. The scalar fields 
are parameterized by 


<f> = 


G+ 

v + £h + iG 0 ) 


y = 


rj-< 


1 


S — -j=(v s + S R + iGs), (II.3) 


-^(Vh + ^a) 

where v is the VEV related with the Fermi constant Gf by v 2 = l/(\/ 2 Gf), and v s is the VEV of 
S which breaks the U( l) M _ r symmetry. In Eq. (III. 3D . G ± , G° and Gs are the Nambu-Goldstone 
bosons which are absorbed by the longitudinal component of the W±, Z and an extra neutral 
gauge boson Z' associated with the I/(1) A1 _ T symmetry, respectively. 

The tadpole conditions for ip H and S R are respectively given by 

dV 


dtp H 

dV 


— v I + 


^Ar + ^As* ) = 0, 


dS H 


~ V S ( Vs + ~R~^S<S> + Vs^S — 0 


(II.4) 


Using the above two equations, we can eliminate and fi 2 s . There is no tadpole condition for y Hl 
because the VEV of inert doublet field 7 / is zero due to the unbroken Z 2 symmetry. 
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The Z 2 -odd component scalar fields, 77 , r]A and t]h , do not mix with the other fields, and their 
squared masses are simply given by 


9 9 x ^ a 

m rj± = + -yAfiV, + — A3, 

U 2 V 2 

m riA = lA + ~2^ Sr l + ^4 — A5), 

v“q r iP‘ 

m Ti H = ^r) + ~2^ Sr l + ~(^3 + A4 + A5). 


(II-5) 

(II.6) 

(H-7) 


For the ^-even sector, two CP-even scalar states (p H and Sh are mixed with each other. Their 
mass matrix, M 2 H , in the basis of (<p H , Sh) is given by 


M 2 h = 


v 2 \\ vv s \ 5 $ 


y vv s X 5 $ 2vgX s 

The mass eigenstates for the CP-even states are given by introducing the mixing angle a by 


(n.8) 



cos a — sin a 
sin a cos a 



(II.9) 


In terms of the matrix element expressed in Eq. (III.81) . the mass eigenvalues are 

m 2 h = cos 2 a(M 2 H )n + sin 2 a{M 2 H )22 + sin 2a(M 2 H ) 12 , 
m 2 H = sin 2 a{M 2 H )n + cos 2 a(M 2 H )22 — sin 2a(M 2 H )i2, 


(II. 10) 
( 11 . 11 ) 


and the mixing angle is 


tan 2a = 


2(M 


h) 12 


( 11 . 12 ) 


(M 2 H )n - (M 2 h W 

We define h as the SM-like Higgs boson with the mass of 126 GeV. Thus, H corresponds to an 
additional singlet-like Higgs boson. Finally, if the conditions 


Ai > 0 , A 2 > 0 , A 5 > 0 , 

(11.13) 

1 1 

As$ + V 2 V^As > 0 , Xs v + —j^V^Xs > 0 , 

(11.14) 

A 3 + ^ VAiAs + min(0, A 4 ± A 5 ) > 0. 

(11.15) 


are satisfied, the Higgs potential Eq. (lII.2| l is bounded from below. 
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FIG. 1: Feynman diagram for neutrino masses at the one-loop level. In the internal fermion line, Nk denotes 
the mass eigenstate of the right-handed neutrinos. 


III. LEPTON MASS MATRIX 


The mass matrices for the charged-leptons and right-handed neutrinos are defined as 


-£mass = (e, n,T)Me(e, /r, r) J 


+ N R> N R)MN{N e R , N£, N^ r ) t + h.c., 


(III.l) 


where e, n and r are, respectively, (e^ + e^), (/x L +/r R ) and (t l + t r ). After the phase redefinition 
of the fields, and N R , the mass matrices can be written in the form 


M t = -^=diag(|y e |, |y M |, |y T |), M N = 


where 0 R is the remaining unremovable phase. Notice here that the U( 1)^_ T symmetry predicts the 
diagonal form of the mass matrix for the charged leptons. The mass matrix Ad tv is diagonalized 
by introducing a unitary matrix V satisfying 


( \M ee \ 

^s_\h 1 

v/^e/d 

ls_\h 1 \ 

V2 |ner| 


^1^1 

0 

\M llT \e ieR 

(HI-2) 

V^l^l 

\M„ T \e ie * 

0 ) 



V t M n V = Ad^ ag = diag(Mi,M 2 ,M 3 ). 

The mass matrix for the left-handed Majorana neutrinos is then calculated to be 


(HI-3) 


(•A A. v )ij — ^ 2 | 2 

fc=l-3 V k 


m m ln 


2 2 
mt, A mt,. 

VA In ■ VA 


~ m m M k M k 


(HI-4) 


If we assume m (j = (m 2 H + m 2 A )/2 3> M|, the neutrino mass matrix can be simplified to be 

(M V )y 


Jv ~ 327 t 2 m 2 (fiVik)Mk(fjVjk) 


0 fc= 1-3 

1 A 5 u 2 


327r 2 m\ 


0 u— 




k= 1-3 


(III.5) 















More explicitly, M u can be written as 
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/ 


My = 


fiM u 

feLMn 


fe fi j. M\ 2 f e frM 13 


\fefrM 13 Uf T M 23 e i9 K 


UfrM 23 e 
0 


(III.6) 




where we reparametrized dimension-full real parameters M,, defined as 


'U o i/v 

Mu = M ee , M 12 = M \ 3 = —f=h eT , M 23 = M mt , 


vs_ 

V2‘ 


(IH-7) 


in the unite of — \ 3 v 2 /(32 tt 2 m^). The structure of matrix, Eq. (IIII.6|) . implies that the U( 1) M _ T 
symmetry predicts the so-called two-zero texture form of the Majorana neutrino mass matrix. 
Fifteen patterns of the two-zero texture form have been discussed in Ref. |29], and our form 


corresponds to one termed “Pattern C”. Because of the two zero texture form, nine neutrino 
parameters, three mass eigenvalues, three mixing angles and three (one Dirac and two Majorana) 
CP-phases, are predicted from five input parameters. In the following, we’ll discuss how we can 
determine all the neutrino parameters by five experimental inputs. 

First, we introduce the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix Dpmns H to di¬ 
agonalize the neutrino mass matrix: 


My = Dpmns diag(mi, m 2l m 3 ) Up MNS , 


(HI-8) 


where mi, m 2 and m 3 are the neutrino mass eigenvalues. The PMNS matrix is expressed as the 
product of two unitary matrices 


T7pmns = UP , 


(III.9) 


where 


U = 


1 0 0 


c 13 0 s 13 e iS 


c 12 s 12 0 

0 c 23 s 23 


0 1 0 


— s 12 c 12 0 

_ 0 ~ s 23 c 23 _ 


_- s 13 e ~ l5 0 c 13 


0 0 1 


P = diag(e ip ,e i,T ,l), (III.10) 


with Sij = sin 6jj and c*j = cos Qij. In Eq. (1III.10D . 5 is the Dirac phase, and p and a are the 
Majorana phases. Using the matrix U, Eq. flIII.8|) is rewritten by 


M u = U diag(mi, m 2 , m 3 ) U T , 


(in.li) 


where fh 3 = m 3 e 2,ip , m 2 = m 2 e Zia and m 3 = m 3 


2icr 
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Second, we obtain the following two equations from the two-zero texture form 
[U diag(mi, m 2 , m 3 ) U 7 ] 22 = [U diag(mi, m 2 , m 3 ) U 1 ] 33 = 0. 

This gives 29] 


(III.12) 


m\ 

m 3 

m 2 

m 3 


c 12 c 13 


c 12 ( c 23 s 23) e * <5 2 s 12 S 23 S 2 3 C 


23 


2 i <5 


S 13 2s 12 C 12 S 23 C 23 (e 2 *‘ 5 + sf 3 ) S 13 (c 32 Si2)( c 23 ' S 23) ei<5 

s 12( c 23 — s 23) el<5 — ^ S 12 S 23 S 23 C 2 


s 12 c 13 


-23 


iS 


(III.13) 


S 13 ^' S 12 C 12' S 23 C 23( e2i<5 + S 13) S 13( C 12 S 12)( C 23 •‘’2 3 ) e * <5 

The ratios of neutrino mass eigenvalues and the Majorana phases are obtained from Eq. (IIII.131) 
as 

mi mi m 2 m 2 1 [mil 1 rm 2 l /TTT i ^ 

Ri 3 = — = — , «23 = — = — , P= xarg — , a = -arg — . (III.14) 

m 3 m 3 m 3 m 3 2 L m 3j 2 |_ m 3 

Using 0 < < 7 t/2 (ij = 12, 13, and 23) and 0i 3 -C 1, we obtain approximate formulae for i?i 3 

and i ?23 as 


R13 ~ 


7^23 — 


1 - 


1 + 


2 cot 012 

cot 2023 COS 5 + ( 

/ cot 012 

sin 0 i3 

V sin 01 3 

2 tan 0 i2 

cot 2023 COS 5 + 1 

(tan 0 i 2 

sin0i 3 

V sin0i 3 


cot 26*23 


cot 2023 


1/2 


1 1/2 


(III.15) 


In order to guarantee m 2 > mi ( i.e., i ?23 > R13), we require cot 2023 cos <5 > 0. In that case, we 
obtain i ? 2 3 > 1, which shows that only the inverted mass hierarchy (m 2 > mi > m 3 ) is allowed in 


our model as already mentioned in Ref. 


2M- 


Finally, we define the ratio of two squared mass difference: 


7-> _ A "4 

rii, = 


|Am|i| 


2 2 

777-2 — m l 

1 2 97 

1777-3 — 777-^| 


From Eq. (IIII. 141) . it can be rewritten in the inverted mass hierarchy as 

2 cot 2012 cot 2023 ~ sin 0i 3 cos 5 


Rv = 


p2 p2 

7 t 2.3 

^ 3-1 


cos 2 0 i 2 2 sin 0i 3 cos d — cot 0i 2 cot 20 23 


We can obtain three mass eigenvalues in terms of Am| l5 R\ 3 and R 23 as 


m 3 = 


V'Ami 


21 


V^23-^ 


mi = m 3 i? 3 i, m 2 = m 3 R 2 3- 


(III.16) 


(III.17) 


(III.18) 


13 


Now, we are ready to determine all the neutrino parameters by using five experimental inputs. 
The best fit (3cr range) values in the inverted mass hierarchy are given as follows [311 ]: 

si 2 = 0.323 (0.278-0.375), s 2 23 = 0.573 (0.403-0.640), s 2 13 = 0.0229 (0.0193-0.0265), 

Am 21 = 7.60 (7.11-8.18) x 10" 5 eV 2 , lAm^l = 2.38 (2.20-2.54) x 10~ 3 eV 2 , (III.19) 
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In the following, we present our predictions using the three sets of input parameters; namely, using 
the best fit values (BF), using the upper limit of the 3a range (±3a) and using the lower limit of 
the 3(7 range (—3a). From two squared mass differences, we can obtain the numerical value 

R v = 0.0319 (BF), 0.0322 (±3a), 0.0323 (-3a), (III.20) 

from Eq. (1III. 161) . We can see that the analytic formula of R v in Eq. (IIII.17D is a function of 6. 
From Eq. (IIII.17D and Eq. (ML 201) . we obtain the Dirac phase 

5 = ±1.96 (BF), ± 2.07 (±3a), ± 0.774 (-3a). (III.21) 


All the negative (positive) solutions for 6 are allowed (excluded) by the experimental data at 95% 


CL 


311 ]. so that we choose the negative solution. We then obtain the ratios as 


m i 

m 3 

7711 

lfl 3 

fh\ 

rh 3 


= 1.39 x e 


= 2.03 x e 


1.91i ^2 

77l 3 

1.52 i 

m 3 


= 1.40 x e 


- 2 . 68 * 


= 2.06 x e 


-2.51* 


(BF), 

(±3a), 


= 1.73 x e -1 ' 19 *, ^ = 1.74 x e 2 ' 78i (-3a), 

m 3 


and the mass eigenvalues and Majorana phases from Eqs. (1III. 141) and (IIII.18I) 

( 7711 , 7712 , 7773 ) [eV] = 0.0583, 0.0589, 0.0420, (p, a) = (0.956,-1.34) (BF), 

(7711,7712, m 3 ) [eV] = 0.0533, 0.0540, 0.0262, (p, a) = (0.759,-1.25) (+3a), 

(77ii,77i 2 ,m 3 ) [eV] = 0.0585, 0.0591, 0.0339, (p, a) = (-0.596,1.39) (-3a). 


Using Eq. (jTTTTSj) , we can get the neutrino mass matrix M\ 

( 0.0408 -0.0186 -0.0378 

M u ~ 


Mu — 


-0.0186 

-0.0186 0 -0.0420 - 0.00631i 

0.0378 -0.0420 -0.006317 0 

/ 0.0249 -0.0228 -0.0410 N 

-0.0228 0 -0.0270 - 0.00352i 

0.0410 —0.0270 — 0.003527 0 

/0.0321 


eV (BF), 


eV (+3a), 


Mu - 


0.0399 


0.0271 


0.0399 0 -0.0344 + 0.002527 

\0.0271 -0.0344 ± 0.002527 0 ) 


eV (-3a), 


(III.22) 


(III.23) 


(III.24) 
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where we performed a phase redefinition so that the phase appears in the (2,3)—component as 
in Eq. (1III.6I) . We thus determine our model parameters by comparing each element of the above 
matrix with corresponding one given in Eq. (1III.6D . 

We note in passing that the lightest right-handed neutrino can be a DM candidate in our 
scenario discussed in this section. The phenomenology of fermionic DM is quite similar to that in 
the Ma’s model [4], and its detailed discussions have been presented in Ref. [5]. 


IV. MUON ANOMALOUS MAGNETIC MOMENT AND LEPTON FLAVOUR 

VIOLATION 

The muon anomalous magnetic moment, so-called the muon g — 2, has been measured at 


32] 


Brookhaven National Laboratory. The current average of the experimental results is given by 

a e * p = 11659208.0(6.3) x 10~ 10 . (IV.l) 

It has been known that there is a discrepancy from the SM prediction by 3.2<r 


Aa M = a® xp - a* M = (29.0 ± 9.0 to 33.5 ± 8.2) x 10 


-10 


33] to 4.1fj Q]: 

(IV.2) 


In our model, the dominant contribution to the muon g — 2 is obtained through the one loop 
diagram where the muon and the extra neutral gauge boson Z' of the U( l)^_ r symmetry are 
running in the loop. The resulting form is given by 

2rx(l — x ) 2 


A a. 




/ dx ——-rg—-—, (IV.3) 

87 U J o r(l — x) z + x 

where gz’ and mz> are the U( l) /i - r gauge coupling constant, the mass of Z', respectively, and 
r = (rn^/mz') 2 . 

On the other hand, the parameter space on mz' and gz' has been severely constrained 
by the neutrino trident producti on p rocess 36] observed in neutrino beam experiments at the 


CHARMII 131 and at the CCFR 


, whose measured cross section well agrees with the SM pre¬ 


diction. For example, gz' > 0.1, gz> > 0.02, gz' > 0.002 and gz' > 0.001 have been excluded with 
95% CL in the cases of mz 1 = 100, 10 , 1 and 0.1 GeV, respectively [36|]. However, we note that 
the muon (g — 2) in our model is not constrained by the dark photon search experiment at BaBar 


because Z' does not couple to the electron in our case 


3 , 


4C] 


By taking into account the constraint from the neutrino trident production, the discrepancy in 
the muon g — 2 can be compensated to be less than 2cr by mz> — 200 MeV with gz' — 10~ 3 1 


1 In addition to the Z' loop contribution, there is a negative contribution to the muon g — 2 from the rf 1 and Ni 
















11 


The lepton flavor violation also arises through the loop in our model. The most stringent 
constraint is imposed by the MEG experiment: B(fi —> ey) < 5.7 x 1CP 13 3f|. The branching 
fraction is written by 

ery) ~ (900 GeV 2 ) 2 x £ A^ Vli V 2 *G ( . (IV.4) 

i=l-3 V m v ± / 

If we take X(*=i -.3 /e//iTj 9 ;U 2 i < 0(1O~ 3 ) with m ± = 0(1) TeV, we can avoid this constraint. 
Therefore, the anomaly in the muon g — 2 can be well explained in the favored parameter region 
suggested from neutrino data and lepton flavour violation data. 

We note that our Z' boson does not couple to the SM quarks, because it appears from the 
U(V)h-t g au g e symmetry; i.e., it has a quark phobic nature. Therefore, any constraints for the Z' 
boson using hadron events such as dijet searches cannot be applied. 


V. CONCLUSIONS AND DISCUSSIONS 


We have constructed a one-loop induced radiative neutrino mass model in the gauge symmetry 
SU(2)l x U(l)y x U(l)fj,- T with the unbroken discrete Z 2 symmetry. In our model, three right- 
handed neutrinos are introduced in addition to the SM, and the scalar sector is composed of two 
isospin doublets, one inert and one active, and a U (1) M _ T charged singlet. 

We have shown that the U (1 )^_ T symmetry predicts a characteristic structure of the lepton mass 
matrices. First, the mass matrix of charged leptons is diagonal in the interaction basis. Second, 
the mass matrix of left-handed neutrinos is in the two-zero texture form if inert scalar bosons are 
much heavier than the right-handed neutrinos. The two-zero texture form of the neutrino mass 


matrix has been intensively studied in Ref. 


29l |. and our model provides a texture with vanishing 


(2,2) and (3,3) elements, corresponding to “Pattern C" in 


29( |. In this pattern, only the inverted 


mass hierarchy is allowed. And we only need five input experimental data to fix the neutrino mass 
matrix. We can choose the most accurately measured ones: two squared mass differences and 
three mixing angles. Using the best fit values of five observables, we obtained non-zero Dirac and 
Majorana CP-phases, and non-zero three neutrino mass eigenvalues. 

We showed that the Z '-loop contribution to the muon g — 2 can explain the discrepancy between 
the current experimental data and the SM prediction if the Z' mass is of 0(100) MeV and the 
U (1) M _ T gauge coupling of 0(1O~ 3 ), which has not been excluded by the neutrino trident production 


(i =1-3) loop diagram. However, it can be neglected due to the assumption -C m;± that provides two-zero 
texture in the neutrino sector. 
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process. The constraint from lepton flavour violation such as g 
parameter space favored by the neutrino data and the muon g — 2. 


ey can be avoided in the 


important features is that the SM-like Higgs boson h which was discovered at the LHC 


O' 

the 

41, 

42] 


can mix with the U{ l)^_ r Higgs boson H via the mixing angle a. As a consequence, the coupling 
constants of h with the SM gauge bosons hVV (V = W, Z ) and fermions hff can be universally 
deviated from those of the SM predictions by the factor cos a. Since the pattern of the deviation 
in the h couplings strongly depends on the structure of the Higgs sector as discussed in Ref. {tJ], 
we can indirectly probe the model by looking at the deviation even if we cannot discover new 
particles such as H. In future collider experiments such as the LHC Run-II, the high luminosity 
LHC and the International Linear Collider (ILC), the Higgs boson couplings are expected to be 
measured quite accurately, especially they can be measured. In particular at the ILC with the 
collision energy of 500 GeV and the integrated luminosity of 500 fby 1 2 3 , the hVV and hff (/ = 6, r 
and t) couplings can be measured with about 0.4% and 0(1)% [4J], respectively. Therefore, we 


can test our model by the comparison between the precisely measured Higgs boson coupling and 
the theory predictions. 
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